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ABSTRACT 

We present Monte Carlo simulations of radiative transfer in magnetar at- 
mospheres. We include the effects of vacuum polarization, electron and proton 
scattering, and free-free absorption. Simulations are performed for the atmo- 
sphere model with the magnetic field perpendicular and also tilted with respect 
to the neutron star surface, and we show that the average spectrum does not 
strongly depend on the orientation of the magnetic field. We investigate the re- 
gion of the parameter space where the vacuum absorption-like feature appears 
in the spectrum and we analyze the shape of the proton cyclotron line. Our 
results indicate that the existence of the vacuum polarization feature should be 
a general attribute of soft gamma-ray repeaters burst spectra, provided that the 
energy release takes place at the sufficiently dense region, and the atmosphere 
scaleheight is large enough. We discuss the existence of such a feature in recent 
observational data on these sources. 

Subject headings: magnetic fields — radiative transfer — stars: atmospheres - 
stars: neutron — X-rays: bursts 
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1. INTRODUCTION 

Soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) belong to a 
group of objects called magnetars, since the main source of energy powering their radia- 
tion comes from the superstrong magnetic field. The first observational suggestion of the 
superstrong magnetic field came from analysis of the March 5th burst properties (Mazets 
et al. 1979; Paczynski 1992), and the basic model of magnetar emission was proposed by 
Thompson & Duncan (1995). The existence of superstrong magnetic fields has been inferred 
from the spindown measurements (Kouveliotou et al. 1998). The AXP population exhibits 
a lot of similarities to the SGRs. Their spins and spindown rates are similar, and recently 
AXPs have also shown bursting behavior (Gavriil et al. 2002). For a review of the properties 
of magnetars see Hurley (2000) and Mereghetti (2001). 

In the conditions prevailing in a magnetar atmosphere, the basic contributions to the 
opacity come from the electron scattering, free-free absorption, and the proton scattering. 
The atmospheres should be completely ionized during the bursts, if they are composed mainly 
of hydrogen, yet in the quiescence some non-ionized material can be present. The density 
changes very rapidly since the gravitational acceleration is in excess of 10 14 cm s -1 . For typical 
magnetar fields, the proton cyclotron frequency ftoj p = 6.4£?i 5 keV, where Bi 5 = _B/10 15 G, 
comes right in the X-ray range. Important additional effects arise from the polarization of 
the vacuum. These effects for strong magnetic fields have been discussed by several authors 
(e.g., Tsai & Erber 1975; Ventura et al. 1979; Pavlov et al. 1980, for a review see, e.g, 
Pavlov & Gnedin 1984; Meszaros 1992). 

Modeling the atmospheres of strongly magnetized sources is usually done using the two 
normal modes description. There are basically two methods to solve the radiative transfer 
equation: the finite difference scheme and the Monte Carlo method. The finite difference 
scheme was used to model the radiation spectra of X-ray pulsars (e.g., Meszaros & Nagel 
1985; Bulik et al. 1992), and also applied to the case of the superstrong magnetic fields 
(e.g., Bezchastnov et al. 1996; Ho & Lai 2001; Ozel 2001; Zane et al. 2001). The Monte 
Carlo method was also used in modeling of accreting X-ray pulsar spectra (Wang et al. 
1989; Isenberg et al. 1998) and to the atmospheres with superstrong magnetic fields (Bulik 
& Miller 1997). The difference scheme method allows to construct self-consistent models of 
static atmospheres, yet requires a careful choice of the grids in order to resolve the sharp 
resonances (see, e.g., Ozel 2003). Imposing the radiative equilibrium is easy in this approach. 
Inclusion of the Comptonization effects is, however, much more difficult. In the case when the 
magnetic field is tilted, not perpendicular to the surface, an accurate solution of radiative 
transfer becomes extremely tedious within the difference scheme. The azimuthal angle is 
an additional degree of freedom, and the dimensions of the matrices increase significantly. 
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However, approximate solutions can be obtained in the diffusion approach (e.g., Pavlov 
et al. 1995). The Monte Carlo method allows to treat the Comptonization effects easily 
and permits to resolve sharp resonances. Imposing the radiative equilibrium constraints is 
much more difficult. Past attempts to use the Monte Carlo technique suffered from the 
limited statistics. The method also allows to investigate the tilted field the length of 

calculations does not increase strongly when the rotational symmetry is broken. It has to be 
emphasized that a comparison of the two methods provides a cross-check to independently 
confirm the results. 

The effects of the vacuum polarization in application to the superstrong fields have been 
first analyzed by Bulik & Miller (1997), who solved the radiative transfer using a Monte Carlo 
method. They modeled the atmospheres with the temperatures above 10 keV, typical of the 
SGR bursts. They have shown that vacuum polarization leads to appearance of a broad 
absorption-like feature in the region around ~ 5 — 15keV in the presence of superstrong 
magnetic fields. The vacuum feature in the spectrum appears because of density dependence 
of the vacuum resonance energy. The radiative transfer is dominated by the extraordinary 
mode because of its much lower cross-section in comparison to the ordinary mode. The 
vacuum resonance then leads to creation of an opaque layer for the low-cross-section mode, 
which effectively changes the depth of its photosphere. Bulik & Miller (1997) neglected 
the effects of proton cyclotron resonance. In recent years, the difference scheme solution of 
magnetized atmospheres has been applied to compute the spectra of magnetars in a series 
of papers by Ozel (2001), Ho & Lai (2001, 2003) and Ho et al. (2003). They have mainly 
considered the case of emission in quiescence, with temperatures around 10 6 K. The proton 
cyclotron resonance was included and shown to appear as a rather prominent line. Proton 
cyclotron resonance has also been discussed by Zane et al. (2001) and Ozel (2003). 

In this paper we present a Monte-Carlo-type calculations of the magnetar atmospheres. 
Our main goals are to include and investigate the effects of vacuum polarization the proton 
cyclotron resonance, calculate the polarization of the outgoing radiation, and solve the ra- 
diative transfer in the case of an atmosphere with tilted magnetic field. In §2 we present the 
calculations of normal mode opacities and polarization vectors, and in §3 we describe the 
Monte Carlo procedure. §4 contains the presentation of the results, which are discussed in 
§5. 

2. RADIATIVE OPACITIES AND POLARIZATION 

Transfer of electromagnetic radiation in a magnetized medium can be conveniently de- 
scribed in terms of normal modes. The complex refractive indices of these modes nj are 
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determined from the wave dispersion equation, which for the medium with the electric per- 
mittivity tensor e and the magnetic permeability tensor \i becomes 

kx[/i' 1 (kxE)] + Q 2 eE = 0, (1) 

where E is the wave electric field. Solutions of equation (1) can be presented using notation 
introduced by Gnedin & Pavlov (1974). In the coordinate system with the magnetic field 
along the z-axis, and the wavevector in the yz-plane at an angle 9 to the magnetic field, we 
define: 

ni = 1 + ^ (e xx + t yy cos 2 9 + e zz sin 2 9 - e xz sin 29 (2) 
~Vxx ~ Vyy cos 2 9 - ii'l sin 2 9 + n~l sin 29), 

riL = - (e xx - £ yy cos 2 9 - e zz sin 2 9 + e xz sin 29 (3) 
+Vxx - Vyy cos 2 9 - /i" 1 sin 2 9 + n~l sin 2^), 

nc = ^(e^cos^ + e^sin^). (4) 
Then, the complex refractive indices are 

rtj = ni ± ^n 2 L + n 2 c , (5) 

where j = 1,2 labels the normal modes. The (real) refraction coefficients of the normal modes 
are given by the real part of rij, and the total absorption coefficients and the cross-sections 
are proportional to the imaginary part of nf £j = (2uj/c)lm(nj) and 

Oj = (2u;/cN e )lm(n j ), (6) 

respectively, where N e is the plasma number density. 

For the normal modes description it is useful to introduce a complex variable b, or the 
two real parameters q and p: 

b = q + ip = —. (7) 

n c 

The polarization vectors in the rotating coordinates e± = 2~ 1 l 2 {e x ± ie y ), eo = e z , are then 
given by 

e{ = — =C j e^{K j cos 9 ± 1), (8) 
v2 

e o = CjKjsm9, 

where 



Kj = b[l + (-l)Vl + 6" 2 ] 
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and Cj = (1 + \Kj\ 2 ) is the normalization. 

The polarization of the modes is determined by the complex quantity 

aj = \ aj \e**, (9) 

related to b through 

a 3 = -r 1 t Vl + b~ 2 . (10) 
The modulus of otj determines the ellipticity 

v \ a j\ - 1 (n) 
3 M + l' 

with \Vj\ being equal to the ratio of the minor axis to the major axis of the polarization 
ellipse, and the sign of Vj determining the direction of rotation of the wave's electric field 
vector. The position angle Xj is measured between the major axis of the polarization ellipse 
and the projection of the magnetic field B on the plane perpendicular to the wavevector k. 

The Stokes parameters for a given normal mode can be written as 

I- = Q. = pj^Ij, Uj = p{jlj, Vj = Pylj, (12) 

where Ij is the intensity of the mode, and p^p^^Py are defined by 

_ 2Re(a j ) , _ 2Im(aj) _ Kf - 1 

If the waves propagate independently, the total Stokes parameters are additive. Therefore, 
for a radiation which can be regarded as a superposition of different waves with intensities 
Ik, they are 

i = j2h,Q = J^pq 1 ^ u = J2pu^ v = J2pv^ ( 14 ) 

fc k k k 

In the Monte Carlo approach this allows one to calculate the Stokes parameters of the 
outgoing radiation by summing up the contributions of the individual photons. 



2.1. Normal modes in the system with cold plasma and vacuum 

To find normal mode properties through equations (2-4) and (5) one must determine 
the permittivity tensor e a b and the magnetic permeability tensor fi a i>. In the medium in 
which the plasma density is low, the contribution of the plasma to the dielectric tensor 
comprises a small departure from unity, so that — 1| 1. The magnetic permeability 
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tensor is essentially nj,^ = 5 a b- Similarly, when the magnetic field is small, the vacuum 
contributions |e^ ac ' ) — 1| <C 1 and \(i~£^ vac ^ — 1| <C 1. Therefore, for such a system the 
contributions of the plasma and the vacuum can be added linearly, and the total dielectric 
tensor and total inverse magnetic permeability are 

, _,W i ,H r -l _ -l(mc) /-,r\ 

Because of this, one can rewrite equations (2-4) in the form convenient for further discussion 
of normal mode properties 

ni = 1 + nf + nT c , n L = n? + n v L ac , n c = ng, (16) 

where the contributions to n^, n?l and ng come from inserting effl — 5 a 6 into equations (2-4) 
in place of t a b ■ 

Due to charge symmetry, the vacuum dielectric and magnetic permeability tensors e^ oc ^ 
and /i~ b 1( ' mc ' > are diagonal and real. Therefore, n™ c = (compare eq. 16), and from equations 
(5) and (16) one can write for the refraction indices in vacuum 

nl ac = ( 17 ) 
n v 2 ac = l+n v j ac + n v L ac . (18) 

The refraction indices n\ a 2 c have been calculated for fku <C m e c 2 and arbitrary magnetic 
field strength by Tsai & Erber (1975): 

nT c = l + ^ m (h)sm 2 9 (19) 
n v 2 ac = i + ^L v± {h)sm 2 6, (20) 

47T 

where ctp is the fine-structure constant, h = B C /2B, where B c = m 2 e c z /he is the critical field 
strength, and 

v± (h) = -4/ilnr(l + / i )+4/ i 2 ^(l + /i) + 2/;,ln/ i + 2/;,[ln(27r)-l]-4/ i 2 + ^, (21) 

2 

r)\\(h) = 81nri(l + /i)-4/ilnr(l + /i)--^(l + /i)-2/iln/i-2/i 2 + 2/iln(27r) 

o 

+^ + i-8L, (22) 
In these expressions, ^f(x) = d/dx\nT(x), 

lnr!(x)= f dt\nT(t) + ^[x-l-\n(27i)}, 
Jo 2 
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and L\ = 1/3 + J* dxlnTi(l + x) 0.249. Using equations (17-20) one can finally write 

nr = ^(nr + nT)-l = ^™ 2 %±W+»7||W] (23) 
nr = l(n? c -nin = ^s™ 2 e[v±(h)- m (h)\. (24) 



The cold electron-proton plasma dielectric tensor is 

= f (25) 



s=e,p 



where the sum runs over electrons (e) and protons (p). In this expression, w pe = (AnN e e 2 /m e ) 1 ^ 2 
is the electron plasma frequency, uj pp = w pe (m e /m p ) 1 / 2 is the proton plasma frequency, and 
n* 6 are the corresponding plasma polarization tensors. The general form of n* b , which takes 
into account both the scattering processes and free- free absorption (but see §2.3) is given 
in §A. For the case of pure scattering the polarization tensors are diagonal in the rotating 
coordinates and from equations (A1-A7) we get 

ic = t - . = (») 

where a = — 1,0, +1, and 7 rs is the radiative width: 7 re = (2/3)(e 2 /m e c 3 )cj 2 for electrons 
and 7 rp = •y re (m e /m p ) for protons, and u ts = u + i7 rs . Since in equation (26) the cyclotron 
frequency is defined as ujb s = (q s B/m s c), the electron cyclotron resonance occurs in U^ a for 
a — +1, and the proton cyclotron resonance occurs in W aa for a — — 1. 

Inserting equations (25-26) into (2-4) and using equation (16), we obtain 



iP l _ \ ^ ^ps 

4 ^ u 

s=e,p 



n\ = 



(1 + C o S 2 6) Uts + sin 2 6— 

u ts - ^Bs W ts 



(27) 



P i 



n L 



s=e,p 



P i cos 9 

\ * ^ps 

9 , , , ,2 , ,2 • ^ Zy ^ 



Damping of particle motion which enters the expressions (A1-A7) is usually very small, 
7||,a) l±,s ^ w ? an d can De neglected in the derivation of normal mode polarization vectors. 
These can be thus obtained from equation (8) by taking expressions (24, 28-29) and (16) 
with 7 rs — > 0, and calculating parameter b in equation (7). 
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2.2. Normal mode polarization properties 

As mentioned above, the normal mode polarization can be described by the complex 
quantity b, which determines the polarization characteristics through oij (eqs 9-11). Calcu- 
lating til and nc in the manner described in §2.1 from the general expressions (A1-A7), and 
neglecting terms proportional to 7jj s and jj_ >s , we obtain 

q « g He^lei^Bp ~ + ^| p (w| e - to 2 ) + u 2 {u 2 Be - iJ 2 )(uJ 2 Bp - u 2 )V (h)] , (30) 

Q0 r 22/2 2 \ 22/2 2\ 

P ~ — [7||,cWpe w Be( w - ^Bp) + 7||, P ^B P (^ ~ W Be ) 

+2c 4 ( 7 x, e (cu 2 - uo 2 Bp ) + 7±, P (cu 2 - u&))V(/i)], (31) 



where we denoted 

sin 2 # 

and 



u) 2 pe u) Be u)(u 2 - u% p ) + ^ 2 p ^ Bp ^ 2 - u% e )] 1 



^) = ^Mfc)-i7ii(fc)]. 

For most of frequencies and angles q 2 ^> p 2 and ctj, Vj and Xj can be approximated as follows 
(Pavlov & Shibanov 1979; Pavlov et al. 1980): 



1 ± y/1 + q 2 

1 



IT 1P 



(32) 



Vj = ± . (33) 

p 7r 1 ± signg . 

* = T 25vTT? + 2 2 S ' g " P ' < 34 > 

Therefore, the modes are nearly orthogonal (xi — 0, \2 — ±7r/2 and Pi = —V2), and when 
|g| <C 1, they are circularly polarized, \V±\ = IP2I — 1, and they are linearly polarized, 
I Pi I = 1 7^*2 1 — for \q\ ^> 1. The latter case is realized for the electron-proton plasma at 
photon energies far from the critical points and in the range w C 0;^ of interest in this 
paper. 

The critical points occur at photon frequencies for which q(uS) = 0. From equation (30) 
one can see that for the system considered here there exist three such points. The first two 
points are located very close to the electron and the proton cyclotron frequencies, and the 
location of the third one, called the vacuum resonance frequency, is approximately 

1/2 



uj v ttuj pe [?7||(/i) -r] ± (h)]J . (35) 
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The vacuum resonance arises when the contributions of plasma and vacuum to the dielectric 
tensor become comparable to each other. Because ujy depends on plasma density, the actual 
frequency at which the resonance occurs varies in the inhomogeneous atmosphere of the 
neutron star. 

The normal mode polarization properties near the critical points are determined by 
the parameter p, which depends on the mechanism of absorption through the corresponding 
damping rates (see eq. 31). Since near these points q 2 <C p 2 , equations (32-34) become 
invalid and must be replaced by other approximate expressions. Thus, for aj we have 



a J = [l±^l-p 2 }/ip, (36) 

and when \p\ < 1 

V 3 = ±^^p Xj = ^signp, (37) 
hence the modes are circular but the polarization ellipses coincide, and when \p\ > 1 

Vj = 0, Xj = ^(signp ± 1) =F \ tan" 1 vV - 1, (38) 

and the modes are linear but the polarization ellipses neither coincide nor are orthogonal. In 
the limiting case where q = and \p\ = 1, the normal modes are completely linear, Vj = 0, 
and entirely coincide, Xi = X2 = =b7r/4, i.e. they are completely nonorthogonal. 

For each of the critical frequencies, the point where the modes are completely nonorthog- 
onal (also called a "mode collapse point") occurs at the critical angle 9 C such that \p(ou c , d c ) \ = 
1. Therefore, as evident from equations (37-38), the behavior of ellipticities and position 
angles near the critical points u c depends on the relationship between 9 C and the photon 
propagation angle 9. At the passage through q = points there is either a change in Xj by 
±7r/2, when 9 < 9 C (\p\ < 1), or the change in the sign of Vj, if 9 > 9 C (\p\ > 1). This means 
that the separation of the modes into an ordinary and an extraordinary one in the presence 
of critical points becomes ambiguous (Pavlov & Shibanov 1979; Pavlov et al. 1980; Bulik & 
Pavlov 1996). 



2.3. Scattering and absorption cross-sections 

The differential cross-section for the scattering on electrons in a hot magnetized plasma 
is given by 

-£^.(u;9 ^ u'9') =rf- J dpf(p)\(e\n\e')\ 2 5(tu + Au-u'). (39) 
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Here r is the classical electron radius, uj, 9, e and a/, 9', e' are the photon energy, the angle 
with respect to the magnetic field B, and polarization vector before and after scattering, 
respectively, and Auj = uj' — uj. The electron momentum distribution f(p) is taken to be 
a one-dimensional Maxwellian for electrons in the lowest Landau level. Averaging equation 
(39) over 0' we obtain 

(^)fM-^')=^ / dpf(p)^\K^\ 2 \e a \ 2 K\ 2 S(u; + Au;-u;'). (40) 

a 

The components of the electron polarization tensor, with first-order relativistic corrections 
in the electron p and photon k, k! momenta are 

n c w = i ^ u\\ 

++ u' + iVBe - pk' + k' 2 /2 - i 7re ' v ; 

n -- = 1 + ■ , (42) 

uj - uj Be - pk - k z /2 - i7 re 

(r) (p + fc/2)(p + fc-fc72) _ (p-U/2)(p-U + k/2) 

00 "~ uj - pk - k 2 /2 - i-f re uj' - pk' + k' 2 /2 - i 7re ' 1 J 

For coherent scattering, i.e. when Auj = 0, these expressions reduce to those given for a cold 
plasma by equation (26). In that case, after integrating (39) over the final states, we obtain 
the scattering cross-sections from one mode to another 

^ = a T J2K\ 2 \K a \ 2 Al (44) 

a 

where 



and ot is the Thomson cross-section. In the transverse mode approximation, the polarization 
modes satisfy the completeness property X^=i2l e ±| 2 = (1 + cos 2 #)/2 and Y2j=i 2 \ e i\ 2 = 
sin 2 9, so that Ylj=i 2 ^± = ^ anc ^ '^2j=i2 J ^i = 1- Using these relations, the total electron 
scattering cross-section of = ofj + cr?- in the cold plasma limit is 

^ = ^5> a | 2 |iCI 2 . (45) 

In the ionized plasma studied in this paper, protons also contribute to the photon 
scattering. However, under condition of thermal equilibrium, proton mean thermal velocities 
are ~ {m v j'm e ) x l 2 times smaller then the electron ones, and for the considered kT = 10 keV 
plasma Comptonization effects, as of equations (39-43), are negligible. Thus, analogously to 
equations (44) and (45), we have for the proton scattering cross-sections 

/ m \ 

tP 



fe) -EKfin^i, (46) 
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and 



As described by equations (44) and (46), a photon can switch between the polarization 
modes in a scattering act. The mode switching probability for the two modes is defined 

erf 



* = ^ (48) 



where again the subscript (s) is for electrons (e) or protons (p). 

The total scattering cross-section is the sum of the electron and proton components, 
namely 

at = ^ + of. (49) 

This expression is essentially equivalent to that given through the optical theorem by equa- 
tion (6), when 7y iS = 7j_ jS = 7 rs (see §A). However, it is important to note that when vacuum 
polarization effects are considered, definition of the modes through the optical theorem be- 
comes ambiguous. Below the vacuum frequency uy, the extraordinary mode is denoted by 
(J = 1) and the ordinary one by (J = 2) in equation (6), whereas this notation is reverted 
above toy. Thus one has to be careful to choose the correct branches of complex roots so that 
the opacities found using the optical theorem (eq. 6) match these calculated directly using 
the polarization vectors. The extraordinary mode is also referred to as the low-cross- sect ion 
mode. 

Free-free absorption is also a significant source of the opacity over the frequencies and 
densities of interest. The cross-section for this process can be found through the optical 
theorem by taking 711 jS = 711 jCS , 7_i_, s = 7_i_, C s i n equations (A1-A7). However, as has been 
recently obtained by Potekhin & Chabrier (2003), these expressions are inaccurate at uo < 
uoBp- To derive the correct expressions for the absorption cross-section one has to take into 
account the effects of the finite proton mass on the process, which allow the absorption 
to occur in proton-proton collisions and also lead to the modifications of the absorption in 
electron-proton collisions due to the proton motion. For kT = 10 keV, the contribution to 
the absorption from proton-proton collisions is negligible and we can write for the free-free 
absorption cross-section 

°i S = a * Yl 7 — , ^7~~7 vTZ — ^\ e a\ 2 9a, (50) 

where the absorption coefficient 

m e L) 4 (7rm e /cT/2) 1 / 2 
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is corrected for stimulated emission, and the dimensionless magnetic Gaunt factors, g± = g± 
and g = g\\, are evaluated according to the formulas given by Meszaros (1992). The effective 
damping in equation (50) is included through 7 a = 7 Q)Ce + 7 re + 7 rp , where the collisional 
damping 7 QiCe = (a a g a /a T )>y re . 

We present the scattering and absorption cross-sections as a function of photon energy 
in Figure la. Vacuum polarization affects the normal mode polarization vectors, which leads 
to redistribution of the total scattering cross-section between the modes: at the vacuum 
resonance frequency Uy, 0\ « a 2 ~ cr T . The absorption cross-section is also enhanced at 
LOy in the low-cross-section mode (for B = and density p = po,i — 220 gem" 3 (see 
below), hjjy ~ 3.38 keV in Fig. la). The vacuum resonance only slightly affects the proton 
scattering cross-section for the density presented. At this density, the proton scattering cross- 
sections for the two polarization modes are of similar magnitude at the vacuum frequency. 
The opacity in the extraordinary mode below lOkeV is dominated by absorption. We note 
that in the wings of the proton resonance {fajJBp — 6.4 keV) the proton scattering cross- 
section dominates over the electron scattering. The electron scattering is resonant at the 
proton cyclotron frequency due to influence of the proton resonance on the normal mode 
polarization vectors. Figure lb shows the mode switching probability as a function of energy 
for the same conditions. The low- to high-cross-section mode switching probability is greatly 
enhanced around the resonances, while far from them it settles at the value of 1/4. 



3. MONTE CARLO RADIATIVE TRANSFER SOLUTION 

To calculate the spectra and polarization properties of radiation emerging from mag- 
netized neutron star atmosphere we use a Monte Carlo method based on Bulik & Miller 
(1997). Assuming an isothermal hydrogen atmosphere with kT = lOkeV and the density 
profile of the form p = poe z / H , we use an adaptive mesh to produce a table of physical 
depth z as a function of optical depth r for the low-cross-section mode, for various values 
of model parameters: magnetic field strength and inclination and the scaleheight H. The 
normal modes polarization vectors are derived in the cold plasma approximation, and the 
scattering and absorption cross-sections are obtained from equations (49) and (50). The 
thermal effects modify the modes for photon frequencies only very close to ujb p , and their 
net result on the scattering cross-sections is to broaden the proton cyclotron resonance (see, 
e.g, Meszaros 1992). Therefore, these effects, although important for the scattering process, 
can be neglected in the optical depth derivations on the grid. The optical depth depends on 
the photon energy and the propagation angle 9, and in the case of a tilted magnetic field on 
the azimuthal angle as well. The changes in the physical depth can be very sharp when 
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approaching the vacuum resonance, and thus the mesh must be carefully chosen. We use 
260 logarithmically spaced zones in photon frequency from hio = 10~ 2M m e c 2 to hio — m e c 2 , 
90 linearly spaced zones in the angle in the range [0, 7r/2] in modeling the case with per- 
pendicular magnetic field (for the model with a tilted field we use 120 linearly spaced zones 
in 9 and <fi in the range [0, ir]), and 276 logarithmically spaced zones in optical depths from 
10~ 2 to 10 4 , which is well sufficient to achieve the accuracy required. 

The transfer of radiation through the neutron star atmosphere is solved with the method 
of Monte Carlo simulations. Photons are injected at a single layer deep in the atmosphere. 
The energies of the photons are selected from a 10-keV blackbody (Planck) distribution, 
and their propagation directions are chosen randomly from a uniform distributions in cos 9 
and 0. At each interaction (scattering or absorption) we generate an optical depth traveled 
by the photon, and from the table of physical depth versus optical depth for given u>, 9 
and 0, we obtain the new physical depth of the photon. If the photon escapes from the 
atmosphere, its energy and propagation direction are kept, together with its polarization 
characteristics, which are calculated in terms of the Stokes parameters (eqs (9-10) and (12- 
13); b in equation (10) being derived at the depth of the last interaction). Otherwise, 
from equations (45), (47) and (50), we calculate the cross-sections and determine which 
interaction process (electron scattering, proton scattering or absorption) occurs. When the 
electron scattering is chosen, we use the differential cross-section of equation (40) to select 
the new energy and the angle of the photon, and determine whether it switches polarization 
modes. In our procedure we do not integrate the cross-sections over the electron momentum 
p in equation (40), but draw the value of p from the thermal distribution. For the proton 
scattering the new photon parameters are selected using equation (46). Finally, in the case 
of absorption, we draw the new energy of the photon from the thermal distribution and 
select in random the new propagation angles. Then the process is repeated until the photon 
escapes from the atmosphere. We impose a reflection-like condition at the bottom of the 
atmosphere: a photon that wanders below the production depth is lost and a new one is 
drawn from a thermal distribution in its place. In a typical simulation we transfer between 
10 9 and 10 10 photons, a significant increase in comparison with 10 5 in the simulations of 
Bulik & Miller (1997). 

In the present simulations, the effects of mode switching are accounted for in the same 
manner as in the work of Bulik & Miller (1997) — the photon is absorbed when switching to 
the high-cross-section mode. For higher photon energies, electron scattering dominates the 
opacity in the high-cross-section mode and the probability Ph^i of switching to the low-cross- 
section mode is very low (see Figure lb). Thus, the photon may require many scatterings 
until it switches its mode again, during which it may substantially change its energy and is 
likely to be absorbed. For lower photon energies, ftw < 10 keV, and in the range of densities of 
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interest, absorption is the main source of opacity for the high-cross-section mode everywhere, 
except at the proton cyclotron resonance. For scattering with the mode switching in either 
vacuum or proton resonance, the probability of returning to the low-cross-section mode in 
the next scattering is high. However, in the case of electron scattering, the fractional energy 
change in the process is much greater than the widths of the resonances and the photon will 
be usually Comptonized out of the resonance, where Ph^i is low. Thus, the applied procedure 
may not be accurate only for the proton scattering at u ~ u Bp , where the process is already 
approximately described because of the thermal effects neglection in the derivation of the 
optical depths. Likewise in Bulik & Miller (1997), the physical depth of the photon is not 
changed at mode switching, because the mean free path is small in the high-cross-section 
mode and the net distance traveled by the photon in this mode is negligible. 



4. RESULTS 

We present our numerical results for the case of the fully ionized hydrogen atmosphere 
with the magnetic field strength B = 10 15 G. This value is chosen in order to conveniently 
discuss the proton effects on the spectral and polarization characteristics of soft gamma- 
ray repeaters. The proton cyclotron resonance thus occurs at Hlob p = 6.4 keV, and the 
departures from the transversality of the normal modes at tko < 2 keV and at the highest 
plasma densities considered, do not influence the results. 

Since, as shown in Bulik & Miller (1997), the shape of the emergent spectrum highly 
depends on the photon production depth, in the following we compare the results obtained 
for two injection depths (in the units of the scaleheight): zo t i = 10H and z ^ = 13if, 
which translates into the plasma density of po,i — 220 gcm~ 3 and po,2 — 4.4 • 10 3 gcm~ 3 , 
respectively, for the assumed density at the surface po = gem" 3 (corresponding to the 
number density N e ~ 6 • 10 21 cm -3 ). We consider two cases of the magnetic field inclination: 
the field perpendicular to the surface and the field tilted at an angle of 45° to the surface. 

Figure 2 shows the example behavior of optical depth r as a function of physical depth 
z for the low-cross-section mode at various photon energies and the propagation direction 
9 ks 48° (cos# = 0.67), for the atmosphere model with scaleheight H = 10cm, and the 
magnetic field normal to the neutron star surface. In this plot, the vacuum resonance 
manifests itself as a sharp increase in the optical depth occurring at the physical depth 
corresponding to the resonant density, at which scattering and absorption cross-sections are 
greatly enhanced at the given photon frequency (see §2.3). The photon energy corresponding 
to the location of the resonant physical depth at z res = zo,i equals hwy,i = 3.38 keV (Fig. 2a). 
For uj > u Vjl the resonant depth increases, reaching z res = z 0t2 at ftwv,2 = 15.15 keV (Fig. 2c). 
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The peculiar behavior occurs at energies near the proton cyclotron resonance. At physical 
depths for which uj v of equation (35) is approximately equal to u>b p , the plasma contribution 
to the dielectric tensor is much greater than that of vacuum and the vacuum resonance does 
not exist, as can be seen in Figure 2b. Because proton scattering and absorption cross- 
sections are large near u>b p , the optical depths are considerably enhanced in this energy 
range. The resonant depths for photon energies shown in Figure 2d lie far below z ,2- 

4.1. Spectra 

Figure 3 shows photon spectra at the top of the atmosphere for the model with scale- 
height H = 10 cm and the magnetic field perpendicular to the surface (B = 10 15 G). When 
the radiation is produced very deep in the atmosphere [z = £0,2), the absorption-like feature 
resulting from the vacuum polarization effects is distinct and broad, and the upper bound on 
the feature is at f\uo V2 ~ 15keV. At energies below and above those at which the absorption- 
like feature forms, the spectrum shows a significant excess over the blackbody flux. For the 
smaller production depth z = zo,i, the overall spectrum is planckian and the absorption-like 
feature is weakly marked and bound in a narrow energy range near fiwyi ~ 3keV. At ei- 
ther production depth, the proton cyclotron line feature at uj pa uj Bp is clearly visible in the 
spectra. 

The absorption-like features in the spectra owe their existence to the influence of vac- 
uum polarization on the absorption and scattering cross-sections, and in the latter case the 
effect of vacuum on the mode switching probability is most important. Optical depths for 
photons interacting in the resonance are large, thus their propagation through the atmo- 
sphere requires several absorption or scattering acts. Absorption can set the photon energy 
out of the resonance, to the region of small optical depths, thus leading to its escape from the 
atmosphere and creating a relative deficit at the original energy. The net effect of scattering 
is similar — at the resonance the mode switching probability is close to unity (see Fig. lb), 
which results in photon absorption at almost every scattering act (see §3). In consequence, 
the effect of Comptonization in shaping the spectrum is negligible in the energy range where 
the absorption-like feature occurs (cf. Bulik & Miller 1997). 

The upper bounds on the absorption-like features in the spectra shown in Figure 3 
correspond roughly to the frequencies u>v,i and uv j2 , at which resonant physical depths are 
equal to the production depths. As can be seen in Figure 2, at energies above hu>v,i the 
optical depths for photons produced at the smaller production depth zo,i = 10H are small. 
Thus most of them are generated above decoupling layer (i.e. at z < Zd ec ) and escape 
freely from the atmosphere. The location of the decoupling layer, estimated roughly as 
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a mean physical depth of the last interaction, is shown in Figure 4, for the representable 
value of the propagation angle 9 = 48°. For energies at which absorption-like feature forms, 
injection depth is on average above the mean decoupling depth, but the differences are small. 
This and the fact that the feature is shallow in this case results from the optical depths of 
photons which form the vacuum feature being only slightly enhanced compared to the regions 
outside of it, so that their propagation is only weakly influenced by interaction processes. 
The situation is different for photons injected at larger depth z 0t2 = 13H, because in this 
case the optical depths of photons forming the feature are very large. As a result, absorption 
and scattering processes modify the spectrum significantly, reducing the photon flux in the 
feature and redistributing it to the regions of smaller optical depths, thus forming "bumps" 
in the spectrum at energies below and above the vacuum feature. Correspondingly, the mean 
decoupling depths in the feature (Fig. 4b) are also considerably reduced. Note also that 
the high-energy part of the spectrum (ou ^> ujb p ) is steeper then the blackbody because the 
cross-section for scattering, that is dominant in this energy range, increases as uj 2 , which 
shifts the distribution of escaping photons to lower energies (Miller & Bulik 1996, see also 
Figs 2d and 4). 

The absorption line at the proton cyclotron frequency fkOBp = 6.4 keV is prominent 
in the spectra. The line is distinct and broader than the natural or thermal width, as 
discussed in, e.g., Ho & Lai (2001). The strength of the line is greater for the larger photon 
production depth z ,2, which can be traced back to the differences in opacities at z 0) i and z 0j2 
at u w u Bp (Fig. 2b). The decoupling density is also reduced at the proton resonance and 
is lower than the decoupling densities in the absorption-like feature (Fig. 4). This results 
from enhanced absorption and the dominance of proton scattering over Compton electron 
scattering near ub p - 

Effects of protons on the plasma response properties are revealed not only as the ab- 
sorption line feature in the spectrum, but also affect its overall shape (see, e.g., Ho & Lai 
2001, 2003; Ozel 2003; Zane et al. 2001). This is exemplified in Figure 3 by short-dashed 
and dot-dashed lines, which show the spectra for the atmosphere model without protons, for 
the production depth z ,i an d ^0,2, respectively. Note that the upper energy bound on the 
absorption-like feature is not changed when the proton effects are neglected, but the shape 
of the feature differs considerably from the model which includes protons. In particular, 
lower energy bound of the feature moves to lower energies in the model without protons. 
The enhanced flux at lower energies in the model with protons occurs mainly because the 
absorption cross-section in the electron-proton plasma is significantly reduced over the pure 
electron plasma value below uj Bp (see §2.3 and eq. 50). 

Figure 5 shows spectra for the model where vacuum polarization is neglected. For the 
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production depth z ,i the flux is essentially blackbody apart from the proton line, which 
shape is not changed compared to the case with vacuum effects included. For photons 
injected at zo^, the striking difference is visible in the proton absorption line, which is much 
broader and stronger when vacuum polarization is neglected. Note that in this case the 
proton line overlaps with the absorption-like feature in the spectrum which includes vacuum 
effects. The reduction in the proton line strength and width by vacuum polarization in such 
a case was also discussed by Ho & Lai (2003) and Ozel (2003). 

As was noted in Bulik & Miller (1997), the exact shape of the spectra and the vacuum 
feature depends also on other parameters than the injection depth, like the magnetic field 
strength, the scaleheight and the pair fraction. Figures 3 and 6-7 show the effects of changing 
the atmosphere scaleheight H . With increasing H the vacuum feature broadens, the upper 
energy bound on the feature remains roughly constant, whereas the lower bound moves to 
lower energies. This is due to the fact that the typical density gradients are smoother when 
H is large and hence the optical depth through the vacuum resonance increases. The depth 
of the proton cyclotron line also increases with scaleheight, which is especially marked for 
the lower production depth zq^. Note that in this case the vacuum feature and the proton 
line are considerably influenced by the value of H and, in particular, they may be virtually 
nonexistent for low scaleheights (e.g. Fig. 6). 

The spectra presented are essentially independent of viewing direction except for the 
angles near ~ 90°. In this case, the vacuum resonance disappears (see eqs 30-31) which 
results in the lack of the absorption- like feature in the spectrum. Further, the depth of the 
proton line is considerably enhanced and the high-energy tail is reduced, which is due to 
geometric effects. 

Figure 8 shows integrated photon spectra for the model with the magnetic field tilted at 
the angle 45° to the neutron star surface (B = 10 15 G) and the atmosphere scaleheight H = 
10 cm. The overall shape of the spectrum does not significantly differ from the one obtained 
in the model with perpendicular magnetic field, also shown in Figure 8 for comparison. The 
angular dependence of radiation can be conveniently analyzed in a coordinate system (x, y, z) 
in which the z-axis is directed along magnetic field. In this system, the vacuum resonance 
disappears when 9 = 90°, where 9 is the wavevector inclination angle to the magnetic field. 
At a given propagation angle 9 and growing azimuthal angle (the azimuthal angle is set to 
zero at the — y-axis and B lies in the yz-plane) , the depth of the proton line and the vacuum 
feature increases and the lower bound on the feature moves to lower energies. This is due to 
geometric effects, since at a given propagation angle 9 and growing <f>, the physical depth to 
the top of the atmosphere increases. 



-18- 



4.2. Polarization characteristics 

Since the low-cross-section mode makes the main contribution to the photon flux emerg- 
ing from the optically thick neutron star atmosphere, the observed radiation is expected to 
acquire strong linear polarization. As noted in §2.2, the polarization pattern of this ra- 
diation is modified around the critical points in the system of vacuum and proton-electron 
plasma, which influence the frequency and angular dependences of normal mode polarization 
characteristics. 

Figure 9 shows an example behavior of ellipticities Vj and position angles Xj of the two 
normal modes in the vicinity of the critical points in a plasma of density po,i (z = zo,i) ■ The 
labeling of the modes is arbitrary. The behavior of Vj and Xj near the first critical frequency 
fouj Cj i = hwy,i — 3.38 keV is different from the behavior at the second critical frequency 
hu) c> 2 = fajJBp = 6.4 keV, because for the photon propagation direction we chose 9 = 78.5° 
(cos# = 0.2) 9 > 9 C)1 = 48° (cosfl Cjl = 0.67) but 9 < 9 C>2 = 89.95° (cos# Cj2 = 8 • 10" 4 ) 
(see §2.2). At the passage through the first critical frequency the sign of Vj changes, but 
the position angles remain essentially the same. The modes are nearly linearly polarized 
(they become exactly linear at u Ci i) except around u c> i, where they acquire weak elliptical 
polarization in the frequency range where \q\ is small and \p\ is close to unity. This also results 
in the small nonorthogonality of the modes around uj c ^. In the vicinity of the second critical 
frequency the normal mode polarization is strongly elliptical. The sign of Vj remains the 
same but the position angles jump by ±7r/2 and coincide at u Ct 2- Note also the jump by n in 
the position angle of the low-cross-section mode at energy close to uj c ^- It occurs at the point 
where parameter q diverges asymptotically (\q\ ^> 1) and has opposite sign at both sides of 
the asymptote while at the same time p crosses zero (compare eqs 33-34). Note that the 
ellipticity is not influenced in this case (see also, e.g., Pavlov & Shibanov (1979) and Bulik & 
Pavlov (1996) for a discussion of the energy dependence of the polarization characteristics). 
With growing density, 9 C corresponding to u c = ub p is approximately constant, whereas 
the critical angle corresponding to the vacuum resonance changes. However, at densities at 
which uj Bp pa ujy the critical angles disappear. 

Figure 10 shows the variation of the linear polarization fraction Fq =< Q > /I with 
photon frequency at the two values of the propagation angle, for the atmosphere model 
with H = 10 cm, magnetic field perpendicular to the surface, and the production depths at 
z = z ,i (Fig- 10a) and z = z ^ 2 (Fig. 10b), respectively. The radiation flux is essentially 
linearly polarized. Small deviations from linearity appear near the critical frequencies and 
especially for photons propagating nearly along the direction of the magnetic field. For the 
case presented in Figure 10a, the decoupling depth is very close to the production depth zq.i 
(see Fig. 4) at all photon energies, and the location of the critical points is approximately 



-19- 



as depicted in Figure 9. Decoupling depths of photons in the model with the production 
depth located at z — z 0j2 span a wide range of z between ~ z 0jl and z 0;2 , which spreads the 
deviation of polarization from linearity over the frequency range corresponding to the width 
of the vacuum feature. 



5. DISCUSSION 

We have presented a series of Monte Carlo calculations of the atmospheres in the 10 15 G 
magnetic field. We took into account the electron scattering, the proton scattering, and 
the free-free absorption. The temperature of the atmosphere is taken to be constant at 
lOkeV. The typical variations of the temperature with the optical depth in the difference 
scheme calculations are weak, T « r 1 / 4 for r ^> 1 (Ozel 2001; Ho & Lai 2001). Therefore our 
assumption of uniform T does not affect the results significantly. The main effect of including 
the temperature profile of this type would be to increase the flux at the high-energy part of 
the spectrum (see, e.g., for the emission in quiescence Ho & Lai 2001, 2003) We investigate in 
detail the effects of the vacuum resonance and the proton cyclotron resonance in the spectra. 
We find that the critical frequencies in the normal modes disappear at the density when the 
vacuum frequency overlaps the proton cyclotron resonance. 

In our treatment of the radiative transfer we have neglected the effects of mode conver- 
sion (Lai & Ho 2002, 2003). Mode conversion may occur when photon traverses a layer of 
varying density and encounters the vacuum resonance. At energies for which the vacuum 
resonance layer is optically thick in the extraordinary mode, the photon will be scattered or 
absorbed, which in our method is taken into account. Thus the mode conversion effects do 
not alter the results in this energy range. In the range of energies where the vacuum reso- 
nance layer is optically thin, the low-cross-section mode photon may change its polarization 
to the high-cross-section one and be absorbed. This effectively increases the optical depth 
through the resonance layer. Thus the major effect of including the mode conversion effects 
would be to extend the vacuum absorption-like feature further to the low energies. This 
effect has also been noted in Lai & Ho (2002). 

We confirm the existence of a wide vacuum feature in the spectra. Bulik & Miller 
(1997) suggested that the feature appears as a consequence of the Comptonization effects. 
We found that it is formed also in an absorption-dominated atmosphere, which means that 
the formation of the vacuum feature is a general property expected to occur in spectra of 
magnetar atmospheres. The vacuum feature appears if the plasma density at the photon 
production depth is large enough. When the scaleheight is increased, the vacuum feature 
becomes more prominent since the optical depth through the resonant layer increases (see 
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the estimates in Bulik & Miller (1997)). We have found that the approximate boundary 
of the density-scaleheight space region where the vacuum absorption-like feature appears in 
the spectrum can be estimated as p in j > 200(i?/lcm) _0 - 26 gcm~ 3 , where p in j is the density 
at which photons are injected into the atmosphere. Thus the vacuum feature should be 
prominent in magnetar spectra provided photons are produced at sufficiently large optical 
depths and the atmosphere scaleheight is large enough. This second condition might be ful- 
filled for SGR bursts, since during the burst the atmosphere is probably lifted by radiation 
pressure and plasma is confined by the strong magnetic field. We do not expect a significant 
absorption-like feature in the SGR spectra in quiescence. In this case the effective temper- 
ature is T pa 10 6 K and the vacuum feature would appear only as a broad depression in the 
exponentially falling high-energy tail, as shown in Ho & Lai (2003) and Ozel (2001, 2003). 

We find that the proton resonance appears as a narrow line on a wider background. 
This line is due to the proton resonance itself. The broad background line appears because 
the presence of the proton resonance affects the normal mode polarization vectors in a broad 
frequency range around it, which in turn influences the opacities. 

We also calculate a model atmosphere with a magnetic field tilted at an angle of 45° 
to the surface normal. The spectrum in this case is very similar to the one obtained for 
the model with perpendicular magnetic field. Such field configurations are important when 
considering radiation emitted from the entire surface of the star. Our result shows that the 
spectrum in the tilted field case can be well approximated by the the spectra calculated in the 
simpler case of the field perpendicular to the surface. This is an independent confirmation 
of the earlier results obtained in the diffusion approximation for the pulsar-type magnetic 
fields (see, e.g., Pavlov et al. 1995). For the magnetar-type fields, the spectra for the field 
tilted at 90° were shown to be similar to the 0° case in the Monte Carlo simulations by Bulik 
& Miller (1997), and in the diffusion approach by Ho & Lai (2001). 

We present the polarization of the outgoing radiation. It is essentially linear and the 
influence of the critical frequencies around the vacuum and proton resonances is very small. 
For the case of the small production depth z = 2 ,i, the effects of the vacuum resonance 
appear as deviations from the linear polarization near ou v , especially at the angles close 
to the direction of B. If the injection depth is large, z = 2^2, the effects of the vacuum 
resonance are spread around the width of the feature, as photons arrive to the observer from 
regions of different densities. The proton cyclotron resonance causes similar distortions from 
linear polarization at ub p , regardless of the production depth. 

SGR spectra in outbursts have been first analyzed in the hard X-ray range above 15keV 
and the typical optically thin thermal bremsstrahlung (OTTB) temperatures were found to 
be in the range of 20 to 40keV (Aptekar et al. 2001). The first analysis of the SGR spectra 
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extending to 5keV have been performed for the case of SGR 1806-20 by Laros et al. (1986) 
and Fenimore et al. (1994). They found that the low energy spectrum was inconsistent with 
the extension of the OTTB spectra to the lower energies. Strohmayer & Ibrahim (1998) 
analyzed the RXTE observations of SGR 1806-20 and came up with a similar conclusion. 
SGR 1900+14 was observed by HETE-2, and Olive et al. (2003) have shown that the low 
energy spectrum is also qualitatively similar. It was not consistent with an OTTB model, 
and was modeled by the two blackbody spectra. Recently Feroci et al. (2004) have analyzed 
several bursts detected from SGR 1900+14 by Beppo SAX in 2001, and found that the 
OTTB model works fine above 15 keV but overestimates the flux in the lower-energy range. 
They tried several alternative models and found the acceptable fits with the two blackbodies 
with temperatures of ~ 3 keV and 10 keV, and the respective radii of 13 and 2 km. We note 
that such a behavior of the spectrum resembles the presence of a broad vacuum feature. Our 
model spectra fall below the underlying thermal distribution around ~ lOkeV. A double 
peaked spectrum as shown in, e.g., Figure 7 could well be modeled as two blackbodies with 
the temperatures differing by a factor of a few. We suspect that these might be the hints of 
the existence of the vacuum polarization feature in the spectra of SGR bursts. However, in 
order to definitely claim the detection of such a feature one would have to perform the data 
analysis using the appropriate model spectra. 

We acknowledge support from KBN through grants PBZ-KBN-054/P03/2001 (JN, TB) 
and 2P03D 00125 (TB). TB is grateful for the hospitality of the Iowa State University, where 
this work was completed. 
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A. PLASMA POLARIZATION TENSORS 

The general form of the plasma dielectric tensor which includes the contribution 
from both the scattering and free-free absorption processes was given by, e.g, Pavlov & 
Shibanov (1979). Rewriting their expressions in the coordinate system with the magnetic 
field along the z-axis, and the wavevector in the yz-plane at an angle 9 to the magnetic 
field, we obtain for the plasma polarization tensors components (s — e,p for electrons and 
protons) : 

xx (w + i7±,.) 2 -w|.' [ } 

lu(uj + i<y ±tS ) [u + i7j|^ cos 2 0(1 + 2 sin 2 9) - i7_L, s sin 2 6 cos 29] 
yy (c + i7||, s )[(c + i7^) 2 -<] ' ( "" ! 
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The above expressions take into account the influence the strong magnetic field has on the 
effective frequencies of electron-proton collisions, which then depend on photon polarization. 
Thus the longitudinal and transverse damping rates are 7|| )S = 7|| jCS + 7 rs and 7j_ iS = 7_l iCS + 
7 rs , respectively, where the radiative widths 7 rs = (2/3)(e 2 /m s c 3 )uj 2 and the corresponding 
collision frequencies 7|| jCS = ((J a 9\\/vT)lrs an d 7_l, cs = (&a,g±/&T)7rs (see also Pavlov & Panov 
1976), with o a given by equation (51) and magnetic Gaunt factors <7im. When only scattering 
processes are considered, 7|| )S = 7j_ iS = 7 rs , plasma polarization tensors become diagonal in 
the rotating coordinates and take the form as given by equation (26). 
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Fig. 1. — (a) Scattering and absorption cross-sections as a function of energy at the density 
p = p 01 ~ 220 gem" 3 (corresponding to z ,i — lOif) and propagation angle 6> = 48° (cos 9 = 
0.67). The solid lines show the total (electron + proton) scattering cross-sections for both 
the high- and the low-cross-section mode. The individual contributions to the total cross- 
section in the extraordinary mode from electron and proton scattering are shown by dotted 
and dash-dotted lines, respectively. Note that the vacuum resonance at hwy — 3.38 keV only 
slightly affects the proton scattering cross-section. The absorption cross-sections for the two 
polarization modes are shown by the dashed lines, (b) The mode switching probabilities: 
low to high mode (dashed line) and high to low mode (solid line). 
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Fig. 2. — The dependence of the optical depth on the physical depth in the atmosphere with 
H = 10 cm and the magnetic field perpendicular to the surface, for the propagation angle of 
48°. The vertical lines denote the photon injection depths in our model simulations, and the 
numbers next to the curves indicate photon energies in keV. Panel (a) shows the region of 
photon energies around fkvy,i, while (c) illustrates the range around huv,2- The energy range 
around the proton resonance Jtwb p = 6.4 keV, where the vacuum resonance disappears, is 
shown in (b). Panel (d) illustrates the energy region above the resonances. 
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Fig. 3. — Model spectra at the top of the atmosphere for the case of the magnetic field 
perpendicular to the surface and the scaleheight H — 10 cm. The dotted line is the input 
Planck spectrum with kT = 10 keV. The solid line represents the case of z — z 0j2 injection, 
while the long-dashed line corresponds to the case of z = zq,i- The spectra for the model 
without protons are also presented for comparison: dot-dashed line is for z = Zo,2, and 
short-dashed line for z = zo,i. 



-28- 




1 10 100 

E [keV] 



Fig. 4. — The mean decoupling depths for the photon propagation angle of 48° and for the 
two production depths we consider: z = z 0> i (a) and z = z 0t2 (b) (H = 10 cm). Note the 
difference in vertical scales. In the case of the z = zq^. injection, the vacuum and the proton 
resonances blend each other. 
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Fig. 5. — Model spectra for the atmosphere with H = 10 cm and the magnetic field perpen- 
dicular to the surface in the case when vacuum polarization effects are turned off: dot-dashed 
line is for z = zo,2, and short-dashed line for z = zo,i- For reference, we also present the 
spectra from Fig. 3 for the model with proton and vacuum effects included. 
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Fig. 6. — Model spectra for the atmosphere with the magnetic field perpendicular to the 
surface and the smaller scaleheight H = 1 cm. In the case of the z = Zo,i injection (long- 
dashed line), the vacuum feature nearly disappears, while for the z = Zo,2 case (solid line), the 
lower bound of the feature moves to higher energy compared with the model with H = 10 cm. 
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Fig. 7. — Model spectra for the atmosphere with the magnetic field perpendicular to the 
surface and the larger scaleheight H = 100 cm. The slowly changing density gradients 
increase the influence of vacuum effects for both production depths: z = zo,i (long-dashed 
line) and z = z 0;2 (solid line). 
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Fig. 8. — Model spectra for the atmosphere with the magnetic field tilted at an angle of 
45° with respect to the surface normal and the scaleheight H = 10 cm. The z = z 0jl case 
(long-dashed line) is practically indistinguishable from the model with perpendicular field 
(short-dashed line). For the larger injection depth, z = zo,2 (solid line), the vacuum feature 
is slightly shifted to the higher energies, as compared to the case with perpendicular field 
(dot-dashed line). 
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Fig. 9. — The dependence of the normal mode ellipticities and position angles on photon 
frequency near the vacuum resonance (left) and the proton resonance (right). The extraor- 
dinary (low-cross-section) mode is the one with x ~ tt/2. Thus, the labeling of the modes 
as extraordinary and ordinary does not follow the continuous lines. The photon propagation 
angle is 9 = 78.5° and z = z 0:1 . 
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Fig. 10. — The deviations from linear polarization, 1 — \Fq\, for the atmosphere with the 
magnetic field perpendicular to the surface and the standard scaleheight H = 10 cm for the 
injection depth z = zo,i (a) and z = zq^ (b). The solid line corresponds to the propagation 
angle 9 « 8°, close to the direction of the magnetic field, and the dashed line is for 9 46°. 



